
% read calibrated level parameters
tab=csvread('../temp/level_param.csv');

param_alpha(1,:)=tab(1,:);
param_psi_m(1,:)=tab(2,:);
param_alpha(2,:)=tab(3,:);
param_psi_m(2,:)=tab(4,:);
param_psi_f(2,:)=tab(5,:);

for m=1:2
    for e=1:2
        
        eval(['psi_m(ind' int2str(m) int2str(e) ')=param_psi_m(m,e);']);
        eval(['alpha(ind' int2str(m) int2str(e) ')=param_alpha(m,e);']);
        if (m==2)
            eval(['psi_f(ind' int2str(m) int2str(e) ')=param_psi_f(m,e);']);
        else
            eval(['psi_f(ind' int2str(m) int2str(e) ')=0;']);
        end
    end
end


